Geodesic flows in rotating black hole backgrounds 

Anirvan Dasgupta* 
Department of Mechanical Engineering and Centre for Theoretical Studies 
Indian Institute of Technology, Kharagpur 721 302, India 

Hemwati Nandan'^ 
CN \ Department of Physics 

^ ■ Gurukula Kangri Vishwavidyalaya, 

■ Haridwar 24 9 4 04, Uttarakhand, India 

'sj- ■ Sayan Kar^ 

. Department of Physics and Centre for Theoretical Studies 

Q-i Indian Institute of Technology, Kharagpur 721 302, India 

^ '. 

bX> Abstract 

I i We study the kinematics of timelike geodesic congruences, in the spacetime geometry of rotating 

■ black holes in three (the BTZ) and four (the Kerr) dimensions. The evolution (Raychaudhuri) equa- 
ls ■ 

^f) • tions for the expansion, shear and rotation along geodesic flows in such spacetimes are obtained. 

in: 

. For the BTZ case, the equations are solved analytically. The effect of the negative cosmological 

I constant on the evolution of the expansion (9), for congruences with and without an initial rota- 

■ tion (wo) is noted. Subsequently, the evolution equations, in the case of a Kerr black hole in four 
^ ■ 

^ ■ dimensions are written and solved numerically, for some specific geodesies flows. It turns out that, 

a I for the Kerr black hole, there exists a critical value of the initial expansion below (above) which 

we have focusing (defocusing). We delineate the dependencies of the expansion, on the black hole 
angular momentum parameter, a, as well as on ujq. Further, the role of a and loq on the time (affine 
parameter) of approach to a singularity (defocusing/focusing) is studied. While the role of ojq on 
this time of approach is as expected, the effect of a leads to an interesting new result. 
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I. INTRODUCTION 



Geodesic congruences (flows) and their features in various spacetime backgrounds have 
been extensively studied by mathematicians and physicists alike, over many years. In the 
context of gravity and general relativity, studies on the evolution (along a chosen geodesic 
flow) of the kinematical quantities (the isotropic expansion, shear and rotation (henceforth 
referred as ESR)) in speciflc spacetime geometries, have been carried out. The ESR evo- 

n n 

lution, as is well-known, is governed by the Raychaudhuri equations |2[. In particular, 
consequences (geodesic focusing) of the Raychaudhuri equation for the expansion scalar are 
crucial while proving the singularity theorems in general relativity s], Q] . 

Recently, we have initiated a programme where we intend to study geodesic flows in 
various contexts (an example outside the realm of general relativity, is, deformable media) 
by solving the evolution equations for the ESR, in various scenarios. Some of our work is 

n n 

available in |5[-[8|. 

In order to push things further beyond our investigations in [3] , here we intend to look at 
geodesic flows in rotating black hole spacetimes. To start with, we restrict to three dimen- 



sions, where the well-known b 
Zanelli (BTZ) line element 
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ack hole solution (with rotation) is the Banados-Teitelboim- 



The BTZ, is a solution of the three dimensional Einstein 



equations with a negative cosmological constant. Further, in four dimensions, we have the 
the famous Kerr black hole l|, |2| which is also important, astrophysically. Such black holes 
are characterised by their mass and spin (rotation rate) and the study of geodesic flows 
in their gravitational flelds can provide useful information about the behaviour of matter 
around strong gravitational fields. 

In this article, we aim towards understanding the kinematics of geodesic flows in the 
above mentioned spacetimes, in three and four dimensions. Our article is organised as fol- 
lows. We flrst review the BTZ spacetime along with its geodesic structure in Section |TT1 
The Raychaudhuri equations in the background of the BTZ black hole spacetime are then 
written in a freely falling Fermi normal frame. Surprisingly, the Raychaudhuri equations 
in this case are analogous to those in the flat background, with stiffness ^ identifled with 
the modulus of the cosmological constant. The analytical solutions are discussed and the 
kinematics of flows is revealed accordingly in Section UTTl Next, the Kerr black hole space- 
time geometry (which is a vacuum, axisymmetric solution of the Einstein equations) with 
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the associated geodesic equations is reviewed briefly in Section HVl followed by the derivation 
of Raychaudhuri equations in Section |Vl In view of the absence of the analytical solutions 
of the (geodesic and Raychaudhuri) equations for the Kerr case, we have solved them nu- 
merically. Based on the evolution of ESR variables for various sets of initial conditions 
with different values of the other parameters involved, the generic features of the geodesic 
flows are visualised in Section |Vl The results are also compared with those obtained for the 
Schwarzschild geometry. Finally, in Section IVIt we highlight the important results obtained 
and indicate the possibilities for future work. 



II. BTZ BLACK HOLE SPACETIME 



Let us begin by writing down the line element for the BTZ black hole in (2 + 1) dimensions 

io|,U, 

ds"^ = -{N^ - r^Nl) + N-'^dr'^ + dcf^ + 2r'^N^dt d(j), (2.1) 
Here, the squared lapse function and the angular shift Nih are given by. 
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(2.2) 



with — oo < t < oo, < r < oo, and < < 2tt. In equation (12. 2p . M and J are, 
respectively, the mass and angular momentum parameters of the BTZ black hole, while the 
parameter / is the radius of curvature of the spacetime geometry, related to the cosmological 
constant (A) as = — 1/A. The lapse function vanishes for the following two values of r. 



r± = I 




(2.3) 



where r± identifies the horizons of the BTZ black hole, which exist only when the following 
conditions are satisfied, 

M > 0; I J| < Ml. (2.4) 

The inner (r_) and outer (r+) horizons shrink with the increasing value of the cosmological 
constant. One may also note here that both the roots of the lapse function coincide for the 
extremal case \J\ = Ml. The first integral of the geodesic equations corresponding to the 



line element (12. ip are given as follows 

i = 



12 



13|, 



iV2 



(2.5) 



± 



A^2 



(2.6) 



(2.7) 



where E and L are integration constants which represent the energy and angular momentum, 
respectively. Here = {t,r,(f)) satisfies the timelike constraint u'^Ui = —1 for the timelike 
geodesies. The equations fl2.5p - fl2.7l) characterise the motion of the test particles in the BTZ 
black hole spacetime. 



III. EVOLUTION OF ESR IN THE LOCAL FRAME 

A. Fermi normal frame and the Raychaudhuri equations 

We construct an orthonormal Fermi frame with m*, {a = 1,2) where 

9ijeae^f} = Sai3 ; 9iju'el^ = (3.1) 

This is the local inertial frame of a freely falling observer. The basis set e^} is parallely 
transported along the geodesic. 

The tensor Bij = VjUi, defined in the coordinate frame, may be rewritten in the local 
frame of the freely falling observer as Bap = Bije^^e'p. Further, Bap can be decomposed into 
its trace, symmetric traceless and antisymmetric parts in the usual way. 



BaP — -^d + aal3 + OJap, (3.2) 



where. 



-^Odap = , aal3= \ , UJal3= \ \ (3.3) 

In the above, 6 is the expansion scalar, a+ and ctx are the components of trace-free shear 
tensor (7^/3 while u denotes the component of rotation in the antisymmetric tensor Ua/B- The 
evolution equation for Baf3 in the Fermi normal frame as discussed above is then given as. 

Since we are working in three dimensions, we can make use of the following relation between 
the Riemann, Ricci tensors and the Ricci scalar, 

Rikjl = QijRkl + QklRij — dkjRil — guRjk — 2 {9ij9kl — 9x19 jk) R (3.5) 



and also the relation Rij = Agij. Consequently, the equation for Bap becomes 

B^^ + B^^BJ = A6^p. (3.6) 

Using f l3.2p . and f l3.6p . the Raychaudhuri equations for the ESR variables turn out to be, 

9 + + 2{al + al-u^) + 2\A\ = 0, (3.7) 

cr^ + 6(7+ = 0, (3.8) 

dx+^(Tx=0, (3.9) 

uj + euj = 0. (3.10) 

The above first-order, coupled, nonlinear and inhomogeneous ordinary differenial equations 
have to be solved, for given initial conditions, in order to obtain the kinematical evolution 
of the ESR along geodesic congruences. 

A subtle point may be noted here. Given the (static) velocity vector field = (t, <p, r) 
in fl2.5p - fl2.7l) one can easily find the ESR variables as functions of r using the definition 
Bij = VjUi. This leads to, for example, 

= ^^u' = ± ^ ^- ^ ^ 7T7^, (3.11) 



and the vorticity w = 0. Physically, this represents the kinematics of a geodesic congruence 
in which all geodesies have a fixed value of the constants E and L. Therefore, this expression 
cannot accommodate arbitrary initial conditions on the velocity field of a congruence, or on 
the ESR variables. The solutions obtained by integrating the Raychaudhuri equations (13. 7p - 
(13.1 op directly, for given arbitrary initial conditions, are, in this sense more general, and 
thereby, underscores the importance of these equations. 



B. Analytical solutions 

The equations (I3.7l) - (l3.10p have been exactly solved analytically (see j^, 0| for the detailec 
procedure). It is interesting to note from equations f l3.7p -( l3?T0|l here and (2.8)-(2.11) in js] 
that, kinematically, the equations in the BTZ black hole background can be mapped to those 
for flows in a fiat space with isotropic stiffness |A| {k in {^). In order to solve this set of 
equations, we first define a parameter / = + — and obtain the following analytical 
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expressions for the ESR for the three different cases depending on the signature of /. 
Case A: / > 

apcos(ci + y) 

^=2[a sin(ci + f ) + 6] ' ^^'^^^ 

r {C2, C3, C4} 

'"^■"■<-"'^ |.sin(c. + f)+tl - 
where = 16|A| and = 16 + 6^. Here we have used { } brackets to express the 
solutions in a compact way for the components of shear and rotation which have a common 
denominator as in equation f l3.13p . The parameter b and the integration constants (ci, C2, C3 
and C4) involved can be given in terms of the initial values of ESR (i.e., OqjCx^q, a^o and 
Uq) as described below, 

+ ^) . (3.14) 




p2 V ^0 21, 



= tan- {^-^\ . (3.15) 

{c2, C3, C4} = {a+0, 0-xo, c^o} (a sinci + 6). (3.16) 

Case B: / = 

^ = 2v^tan[v^(ci -t)], (3.17) 

{cr+,a-x,w} = {c2,C3,C4}sec^[A/|A|(ci - t)]. (3.18) 
The integration constants in terms of the initial ESR for this case are as given below, 

ci = tan-i ( — ^ ) , (3.19) 



|A| V^vl^ 

{c.,C3,C4} = i^±^^^^. (3.20) 
sec2(^|A| ci) 



Case C: / < 



apsin(ci + y) 
2 [a cos(ci + f ) + 



(3.21) 



{C2,C3,C4} 

[a cos(ci + f ) + 6] 



= ^— — y— (3.22) 
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Here b is the same as in equation fl3.14p for the case / > and the integration constants in 
terms of initial values of the ESR are, 

200 



c\ = tan ^ 



{c2, C3, Ci} = {cr+o, o'xo, uJo}{a cosci + b). 
All the solutions mentioned above become singular in finite time. Note, however, that this 



(3.23) 
(3.24) 





FIG. 1: Effect of A on the evolution of the expansion scalar 9{X) for the cases (i) without initial 
rotation, and (ii) with initial rotation (a;o = 0.2). 



>- 



(ii) 




>- 




-2 2 4 6 



FIG. 2: Orbits around the BTZ black hole without initial rotation for (i) with horizon (A = —0.4), 
and (ii) without horizon (A = —0.45). The dotted and dashed lines indicate, respectively, the outer 
(r+) and inner (r_) horizons in the left panel above. 



is a singularity of the congruence. One may note that the kinematics of geodesic flows, as 
governed by the Raychaudhuri equations, in the freely falling frame, is independent of the 
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FIG. 3: Orbits around the BTZ black hole with initial rotation (wq = 0.2) for (i) with horizon 
(A = -0.4), and (ii) without horizon (A = -0.45). 

value of the parameters M, J, E and L. However, the geodesies are dependent on these 
parameters and we use the following values of parameters M = 1,E = 2,L = 2 and J = 1.5 
for visualising the orbits around the BTZ black hole. For the kinematics, we consider two 
cases with A = —0.40 and A = —0.45 which correspond to the BTZ black hole with and 
without horizon, respectively. For these two cases, we observe the effect of rotation on the 
expansion scalar 6 in FigiH The corresponding orbits are presented in Figl2] (without initial 
rotation) and Figl3] (with initial rotation). It is observed that the geodesic congruence is 
necessarily singular without initial rotation (see FigH]), a fact which can also be inferred 
from the analytical solutions. Geodesic focusing always occurs earlier with the shrinking 
of the horizon of the BTZ black hole, due to an increase in the value of |A|, as shown in 
FigUti)- However beyond a critical value of the initial rotation, the congruence singularity is 
removed irrespective of the presence/absence of the horizon, as shown in Figll](ii). It may be 
noted that the FiglT](i) corresponds to the exact solutions of the ESR variables represented 
by equations f l3.12p and (13.131) for the case / > 0. Further, FigHJ^ii) corresponds to the exact 
solutions of the ESR variables represented by equations (13.211) and (I3.22p for the case / < 0. 

It might seem that the expression for 6 in (13.111) does have a dependence on the pa- 
rameters M, J E, L among others. However, one can show using the solutions of the 
geodesic equations (i.e., r{X)) in (13. lip that the final result coincides with one of the solu- 
tions obtained above in this section. Thus, it is surprising to note that, while the geodesic 
trajectories depend on the parameters M, J, etc., the ESR variables depend only on |A|. 
This is essentially an effect of the spacetime dimension (2 + 1 here) where the Weyl tensor is 
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identically zero and no explicit dependence on M, J appear in the ESR evolution equations. 



IV. THE KERR BLACK HOLE SPACETIME 



In this section, we quickly recall the well-known Kerr black hole spacetime (in the Boyer- 
Lindquist coordinates, which are analogous to the Schwarzschild coordinates for a non- 
rotating black hole) with mass M and angular momentum J as characterized by the following 
stationary and axisymmetric metric, 

ds^ = -[l — dt^ d(l)dt + ^ dr^ + dijj^ 



J p2 ^ ' A 

+ ir^ + a^ + J sm^ d<j? , (4.1) 

Here a = J/M, p^ = r"^ + a'^cos'^ip and A = — 2Mr + a^. With a = 0, the Kerr metric 
(14. ip reduces to the Schwarzschild. The geodesic equations corresponding to the metric 
(14. ip are well-known (we do not repeat them here). However, analytical expressions of first 
integrals of the geodesic equations are not known for the general case (except for the choice 
of = IT /2 14]). In our work here, we are interested in timelike geodesic congruences for 
which the constraint u°'Ua = — 1 is satisfied, i.e.. 



goot" + gnr^ + ^22^' + ^33<^' + 'igosi^ + 1 = 0. (4.2) 

Once the initial conditions on the velocity vector components are specified in such a way 
that this constraint is satisfied, it will be satisfied throughout the geodesic evolution. 



V. ESR EVOLUTION 



A. The Raychaudhuri equations 

Consider a congruence of geodesies with the associated time like vector field u", tangent 
to the geodesic at each point and satisfying equation (14. 2p . A transverse metric ha/s on a 
spacelike hypersurface can be written by decomposing the spacetime metric (gais) as defined 
through equation (14. ip into the longitudinal {—UaUp) and transverse part (i.e., /i^/s), 

hafi = gafi + UaUfj] («, /3 = 0, 1 , 2, 3) . (5.1) 
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The transverse metric satisfies u^ha/B = which means that ha/s is orthogonal to and the 
transverse spacehke hypersurface represents the local rest frame of a freely falling observer 
in the given spacetime. The evolution of the ESR on such a transverse hypersurface, can 
be investigated by using a tensor Bajs, which can be decomposed into its trace, symmetric 
traceless and antisymmetric parts as follows, 

Ba/B = -0 + CTa/S + (5.2) 

where 9 = B"^ is the expansion scalar, aaj3 = B(^ai3) — Oha/^/S the shear tensor, and 
Uai3 = B[ai3], the rotatiou tensor with the brackets () and [] denoting symmetrisation and 
antisymmetrisation, respectively. The shear and rotation tensors also satisfy /i"'^ cTq,^ = 
and h"^ ujap = 0. The evolution equation for B^p takes the following form, 

BaP + Ba-yB'^i^ = —Rar)l35 u'^U^ , (5.3) 

Using the equation (15. 3p and following standard methods [l|, [Zl, [Oj, the Raychaudhuri equa- 
tions for the ESR variables in the Kerr black hole spacetime, given by equation (14.11) turn 

out to be, 

e + ^e' + a^-tu^ = 0, (5.4) 
2 1 

2 , , 

(^a(3 + -0 + <Ta ^7/3 + = 0> (5-6) 

where cx^ = (Jq,/3 0"°^, u"^ = ujapoj'^^ and Caiir,s = Ra-qps is the Weyl tensor for the present 
case. One can notice the absence of the geometric part (i.e., Rapu^u^ ) in equation (15. 4p 
since the Kerr spacetime is Ricci flat. Further, the nonzero Weyl in equation (5.5), with its 
dependence on the rotation parameter a (and also M) will lead to a dependence on a in the 
evolution of the ESR. 



B. Visualisation of ESR evolution 



In the absence of analytical solutions for the ESR variables in this case, we have directly 
solved (15. Sp numerically. In order to represent these kinematical quantities at any point, we 
consider a freely falling (Fermi) normal frame with the basis vectors i?^, /i = 0, . . . , 3 (with 
Eq = u") which are parallely transported. To construct such frames numerically, we solve 

10 



the differential equations = (with initial conditions of an orthonormal frame) 

simultaneously with (15 .Sp The tensor Baj3 in the Fermi basis may be represented as, 

111 

Bap = (-9 + au)ey^ + (-6 + ^22)6^6^ + (-9 - an - a22)elel + (ai2 - ^3)6^6^+ 

((721 + W3)e^eJ + (ai3 + W2)e^e^ + (0-31 - uj2)e\el + ((123 - Wi)e^e^ + {a^2 + Wi)ej?,e^ , (5.7) 

where are co- frame basis which satisfy the relation e(^i?" = 5^'. The ESR variables can 
now be constructed from the evolution tensor (15 .Zp . using the basis vectors, as. 





t' — -Da/3 








0"11 = -Ba/3 -E" 


- -9, 
^ 3 ' 


(5.9) 




(T22 = Bai3 E2 


E^ - -e, 
2 3 , 


(5.10) 




= ^{BapE^ E^ 


+ BapE^ E'^), 


(5.11) 


<7l3 


= -{BapE" E^ 


+ Bap E^ Ef), 


(5.12) 


'^23 


= ^{BapE^ E^ 


+ BapE^ E^), 


(5.13) 


Wi 


= -{BapE^ E^ 


- Bap E2 E^), 


(5.14) 


^2 


= -{BapE^E^ 


- Bap E^ ), 


(5.15) 


W3 


= -{BapE^ E^ 


— Bap E^ E2). 


(5.16) 



In order to understand the focusing/defocusing behaviour of a timelike geodesic con- 
gruence, let us redefine the expansion scalar as ^ = 3F/F where the dot indicates the 
derivative with respect to the affine parameter A. The equation (15. 4p can now be written in 
the following Hill-type equation, 

F + HF = 0, (5.17) 

where if = (cr^— a;^)/3 wither^ = 2(cr^^+cr|2+o"i2+o"i3+o"23+^ii^22) andw^ = 2(a;^+a;|+u;|). 
If F — )■ in finite time, we have a finite time singularity in 9 with focusing (defocusing) 
if F < (F > 0). The signature of the invariant quantity H is useful in understanding 
the individual roles of shear and rotation in the occurrence of geodesic focusing/defocusing. 
When H is positive definite (i.e., > w^), there exists conjugate points and geodesic 
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FIG. 4: Defocusing (for > ^c) for the cases (i) without initial rotation, and (ii) with initial 
rotation (wq = 0.1). The effect of (iii) angular momentum, and (iv) initial rotation on time of 
approach to singularity (As). The initial shear in all cases is zero. 

focusing/defocusing takes place. On the other hand, no finite time singularity exists for an 
initially non-contracting congruence (i.e., 6q > 0) if H is negative definite. For 6q < 0, there 
exists a critical value below which focusing/defocusing will take place. 

In the following, the evolution of the ESR variables in the Kerr black hole background, 
under different conditions, are presented and compared with those in the Schwarzschild 
background for some interesting cases. 

We consider the case when if) = tt/2 (equatorial section). We have observed that for an 
initially expanding congruence, geodesic defocusing takes place (i.e., — )• oo) and there exists 
a critical value of the expansion scalar 9c below which there is a focusing (i.e., 9 — ?■ — oo). 
It is important to mention that with the change in the initial conditions on shear and 
rotation, the critical value of the initial expansion scalar will also change. The results 
with no initial shear and rotation are presented in Figs. Il](i) and (ii) for 9q above the 
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FIG. 5: Focusing (for 6q < 9c) for the cases (i) without initial rotation, and (ii) with initial rotation 
{ojQ = 0.2). The effect of (iii) angular momentum, and (iv) initial rotation on time of approach to 
singularity (Ag). The initial shear is ctq = 0.28. 

critical initial expansion (defocusing), and Figs. [5][i) and (ii) for 9q below the critical initial 
expansion (focusing) for both Kerr and Schwarzschild (a = 0) black holes. In order to see 
the effects of the angular momentum a and the initial rotation ujq of the congruence, we have 
calculated the time to singularity for the super-critical and sub-critical values of 9q and 
presented them, respectively, in Figs. IH^iii) and (iv), and Figs. El^iii) and (iv). For the case 
in Fig. m we observe that both angular momentum and initial rotation of the congruence 
assist defocusing. On the other hand in Fig. |5l it is observed that with increase in wq, 
the congruence focuses at a later instant, as expected, while an opposite trend is observed 
for increase in a. Thus, we conclude that the angular momentum a advances the time to 
singularity (irrespective of its type), while the initial rotation ujq of the congruence advances 
defocusing and delays focusing. For this reason, for the same set of initial conditions, 
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focusing/defocusing occurs earlier in the Kerr black hole as compared to the Schwarzschild. 
We have also noticed that the effect of initial shear on geodesic focusing and defocusing for 
the Kerr case is the same as that for the Schwarzschild (results are not presented here). 
Focusing occurs at an earlier instant of time when we have some initial shear in comparison 
to the case with zero initial shear. 

Since all values of M and a do not correspond to a Kerr black hole and the horizon exists 
only when a < M, the limiting value of the angular momentum is a = M (extreme Kerr). 
We have observed that for the extreme Kerr, all the features of the geodesic focusing and 
defocusing remain qualitatively similar as those for the non-extreme Kerr case, though with 
a change in the scale. 



VI. SUMMARY AND CONCLUSIONS 



In this article, we have investigated the kinematics of of timelike geodesic congruences in 
the background of the BTZ and Kerr black hole spacetimes. The important conclusions are 
summarised below: 

• Apart from the cosmological constant (A), the other parameters of the BTZ black hole 
spacetime, do not have any direct role in the kinematics of the geodesic congruence. 
However, the geodesic trajectories surely depend on the black hole parameters. 

• An increase in A in the BTZ case, which implies a larger negative curvature, enhances 
focusing of the congruence. 

• The notion of focusing remains qualitatively similar in the BTZ black hole background 
for the cases with or without horizons. 

• The angular momentum parameter a of the Kerr black hole assists both focusing and 
defocusing when the respective cases occur. Thus, physically, the larger the spin of 
the black hole, the earlier the focusing/defocusing of the geodesic congruence occurs. 
This effect is absent in the BTZ black hole case. 

• In contrast with the role of a, an increase in the initial rotation, uq, of the congruence 
assists defocusing but delays focusing in both the (2 + 1) and (3 + 1) dimensional 

n Q 

rotating black holes. This observation is similar to results obtained previously [5[-[7J. 
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The tidal distortion due to the presence of mass nearby, or accretion of matter into a 
rotating black hole deserves careful attention. The kinematics of inflows can be conveniently 
studied using the Raychaudhuri equations in these backgrounds. The effect of the black 
hole angular momentum parameter on the focussing/defocussing of geodesic congruences, 
as has been observed in this work, can have an interesting role in the accretion process 
and may produce observable effects. Further, for a more comprehensive understanding, the 
kinematics of flows in maximally extended coordinate systems, for black hole spacetimes, 
would be interesting to investigate. Finally, the study of kinematics of null geodesic flows 
in such spacetimes is surely worth trying out in future. 
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